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Abstract 

Using data collected with the OPAL detector at LEP, we have searched for the pro- 
cesses e + e~ — > Z° — > pe~,p/x~ and the charge conjugate final-states. These would violate 
the conservation of the baryon-number B, lepton-number L and the fermion-number n = 
(B+L). No evidence for such decays has been found, and the 95 % confidence level upper 
limits on the partial widths r(Z° — ► pe) and r(Z° — ► p/x) are found to be 4.6 and 4.4 keV 
respectively. 
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1 Introduction 



No symmetry principle is known in physics which may guarantee the conservation of baryon (B)- 
or lepton (L)-number. Nevertheless, there are two important observations which are relevant to 
B and L non- conservation in a very profound way: namely baryogenesis Jl| , and the realization 
that B and L are not strictly conserved in the Standard Model of electroweak interactions 0. 

In the standard hot big bang model || , one of the most obvious relics from the early universe 
is the baryon. But the universe, at some distance scale around us, £b, is practically 100 % charge 
asymmetric, with baryon number density very much exceeding that of antibaryons, ub ^> n B 
(and n e - 3> n e +). The distance Ib is at least of the size of our visible universe. The standard 
electroweak theory contains in principle all the elements necessary for the generation of the 
baryonic asymmetry of the universe @, |5|. These theoretical elements are 

(i) anomalous electroweak baryon- number non-conservation [0], 

(ii) C- and CP-violation in the fundamental gauge and Higgs interactions of the quarks, 

(iii) deviation from thermal equilibrium, assuming the cosmological electroweak phase transi- 
tion is first order ||, which occurs at temperatures of about ~ 100 GeV.Q 

At low temperatures, this anomalous electroweak B violation is negligible |l|, For in- 
stance, the proton lifetime, determined from the absence of the decay p — > e + Z°* — > e + e + e~ 
was found to be larger than 5 x 10 32 years at 90 % confidence level J7|. On the other hand, 
at the Z° resonance with an equivalent temperature of 90 GeV, the decays Z° — > pe/p/i may 
conceivably occur at an observable level due to the anomalous interactions mentioned above. 

In this paper, we present the results of a search for the decays Z° — > pe~,p/i~ and their 
charge conjugate final-states, pe + and p/x + , with the OPAL detector at LEP. These decays 
would violate B, L and the fermion number, n = (B+L) while conserving (B-L). We have 
searched for events characterized by a pair of back-to-back charged particles, where one is 
identified as a lepton and the other as a proton. Although other rare decay modes of the Z° 
can also be handled in the same analysis, they are not discussed here since they either violate 
conservation laws attributable to some fundamental symmetries, or since stringent constraints 
on their rates have already been set by OPAL and other LEP experiments ||. 

2 The OPAL detector and data sample 

A complete description of the OPAL detector can be found in Ref. 0, and only those com- 
ponents most relevant to the identification of high energy electrons, muons and protons are 
described briefly here. 

1 Sakharov suggested Q that baryogenesis took place immediately after the big bang, at a temperature not 
far below the Planck scale of 10 19 GeV. This scenario was explicitly realized with the advent of the Grand 
Unified Theory (GUT), which predicts B and possibly CP violations, mediated by very heavy gauge bosons, X 
and Y. The GUT baryogenesis is appealingly simple but does not easily fit into an acceptable cosmology. 
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The central detector provides charged particle tracking over 96% of the full solid angle in 
a 0.435 T uniform magnetic field parallel to the beam axis |. It consists of a two-layer silicon 
microstrip vertex detector, a high precision drift chamber, a large volume jet chamber, and 
a set of z chambers measuring the track coordinates along the beam direction. The sensitive 
volume of the jet chamber is a cylinder of about 2 m in radius and 4 m in length. It is 
divided into 24 sectors, each equipped with 159 sense wires and two cathode wire planes. Up 
to 159 measurements of position and charge deposition per track are possible. The momentum 



resolution of the whole tracking system is cr pt /p t ~ y0.02 2 + (0.0015 • p t ) 2 where p t is the 
momentum component transverse to the beam direction. The charge measurements were used 
for the calculation of the specific energy loss dE/dxoi a track. The dE/dx resolutions obtained 
with a number of samplings iV > 130 are = 3.1 % and 3.8 % for muons in Z° — > fi + n~ and 
for minimum ionising pions within jets, respectively fllQ"| . A detailed description of the main 
features of the jet chamber can be found in . 



The magnetic coil is surrounded by a time-of-flight counter array (TOF) which measures the 
time-of-flight for a charged particle from the interaction region with a time resolution of about 
300 ps and can be used for the rejection of cosmic ray muons. A lead-glass electromagnetic 
calorimeter (ECAL), located outside the magnet coil, covers the full azimuthal range with 
excellent hermeticity in the polar angle range of | cos#| < 0.984. In the barrel region (| cos#| < 
0.82), the ECAL is composed of 9440 lead-glass blocks, each 37 cm in depth. The blocks 
are approximately 10 x 10 cm 2 in cross section, and the calorimeter is typically 24.6 radiation 
lengths deep. The energy resolution of the barrel ECAL is approximately <Je/E ~ 2.3% 
for 45 GeV electrons. The magnet return yoke is segmented and instrumented for hadron 
calorimetry (HCAL) and consists of barrel and endcap sections along with pole tip detectors 
that together cover the region |cos#| < 0.99. The HCAL provides tracking for muons above 
3 GeV and sampling for the showers produced by hadrons interacting with material (a total of 
~ 8 absorption lengths to the back of HCAL in the barrel region) inside the lead glass of ECAL 
and the return yoke of HCAL. Four layers of muon chambers cover the outside of the hadron 
calorimeter. 

The electromagnetic calorimeters close to the beam axis complete the geometrical accep- 
tance down to 24 mrad in 8, except for the regions where a tungsten shield is present to protect 
the detectors from synchrotron radiation. These are lead-scintillator sandwich calorimeters and, 



at smaller angles, silicon tungsten calorimeters JL^] located on both sides of the interaction point 
and used for luminosity measurement. 

The data used for this analysis were recorded at LEP during 1990 - 1994. The integrated 
luminosity of the data sample is 119 pb _1 [0| corresponding to 5.2 million Z° decays. 



3 Monte Carlo simulations 

The simulation of Z° — > pe~,p/i _ decays and their charge conjugate final-states begins with 
20,000 Z° — > events that were generated at the Z° resonance peak with ^Js = 91.2 GeV 



using the KORALZ WM Monte Carlo. The masses, particle types and the four momenta vec- 



tors of the muon pairs were then changed into those of pe , pe , p/z or p/i + pairs with equal 

2 The OPAL coordinate system is defined so that the z axis is in the direction of the electron beam, the 
x axis is horizontal and points towards the centre of the LEP ring, and 9 and <f> are the polar and azimuthal 
angles, defined relative to the +z- and +£-axes, respectively. The radial coordinate is denoted as r. 



5 



probability and the necessary corrections were made to preserve momentum and energy con- 
servation. Responses of the OPAL detector to these particles were then obtained by processing 
the four vectors through the OPAL Monte Carlo simulation chain ||15|| . 

The analysis procedure consists of the selection of a pair of collinear charged particles in 
the barrel region of the jet chamber, each with approximately the beam energy, followed by 
particle identification. Major backgrounds were anticipated from mis-identification of the high 
energy charged leptons in Z° — > e + e~, and r + r~ decays. Large samples of Standard 

Model lepton pairs ( 600 000 942 000 e + e~ and 447500 r + r~ ) were generated with the 

KORALZ and BABAMC [16 Monte Carlo programs. In addition, 1.3 million Standard Model 
Z° decays into qq's, where q = u,d,c, s and b quarks, were also generated with the JETSET [|17 
Monte Carlo program. All events are then processed through the OPAL analysis chain [15]. 



4 Preselection 

In this analysis, candidate events were chosen from a sample of low (charge) multiplicity events. 
We demanded that the events contained two, and only two, well-measured high energy tracks, 
satisfying the following preselection criteria to guarantee a high degree of redundancy in trig- 
gering and to give almost an 100 % trigger efficiency for the detection of the signal ]TB|, |T9"] : 

• Well-measured tracks were defined by: 

The distance of closest approach to the beam axis in the plane perpendicular to the beam 
axis must be less than 2 cm. 

The distance to the interaction point along the beam axis, at the point of closest approach 
in the plane transverse to the beam axis, must be less than 40 cm. 

The x 2 /n.d.f. of the track fit in the r — <fi plane must be less than 2. 

The first sense wire in the jet chamber used for track fitting must lie within a radius of 
75 cm from the colliding beam axis. 

There must be at least 100 hits used for the ionization energy loss, dE/dx, measurement 
and 4 hits in the z chambers. 

The track polar angle must satisfy | cos#| < 0.7. 

The track must have either its measured p t > 15GeV/c, or its associated Eecal > 
35 GeV. 

• To select the required event topology, the two tracks were required to be: 

oppositely charged and collinear within 59 aco i < 10°, 

in coincidence at the TOF counters, where we demanded that: 

\tisin6\ — t2sin6 2 \ < 10 nsec 

where £x,2 an d #1,2 are the measured time-of-flights and polar angles of the two preselected 
particles respectively. 
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At this stage 139 328 events passed the event preselection criteria. Possible contamination 
from two photon interactions were efficiently removed by the high p t and E EC al cuts. Cosmic 
muons passing through the interaction region were rejected by the coincidence requirement on 
their TOF measurements. The surviving events, consistent with the number expected from 
standard model Z° decays to e + e~, and t + t~, constituted the major background to the 

signal. Comparisons of the data with these Standard Model Monte Carlo simulated Z° decay 
backgrounds are shown in row 3 of Table [I] for the numbers of events that survived preselection, 
and in Figures 1 and 2 for the E tot = E EC al + E HC al and dE/dx distributions of the charged 
tracks, respectively. The Monte Carlo simulated background reproduced the data very well. 

5 Particle identification 

We demanded, among the two pre-selected high energy particles in each event, that there be 
one identified charged lepton and one charged track identified as a proton and not as a lepton. 
Particle identification was achieved using five independent quantities, namely momentum (p), 
ionization energy loss (dE/dx), energy deposits in the calorimeters, Eecal and Ehcal, and 
the information on matching of muon chamber hits to the charged track. We defined a 'tracking 
road' along the direction of each charged track in the OPAL calorimeters with a width: 

20 mrad + 100 mrad/p ( GeV/c), in HCAL 

100 mrad + 100 mrad/p ( GeV/ c), in the muon chambers. 

• Electron identification: 

Electron candidates were defined as those tracks having Eecal > 35 GeV with Eecal /v > 
0.8, and no hits found in the 'tracking road', except for the first layer of the HCAL. 

• Muon identification: 

As muon candidates we considered those tracks which failed the electron identification. 
Those tracks must have its measured momentum, p, satisfying 40 GeV/c < p < 60 GeV/c, 
and have more than 4 hits along the 'tracking road'. To distinguish muons from protons, 
we further demanded that either at least two of the hits on the 'tracking road' be in the 
muon chamber, or there be at least 4 of the hits in the HCAL with no more than one hit 
per layer, or that there be at least 2 of the hits in the 3 outermost layers of the HCAL. 

• Proton identification: 

We considered those tracks as proton candidates which failed the electron identification. 
The proton track must have 40GeV/c < p < 60GeV/c and have at least one hit found 
in HCAL with the expected total calorimeter energy associated with the measured track 
E to t > 10 GeV, but requiring that Eecal < 35 GeV. 

• Event identification: 

Finally, events identified as e + e~ or pairs were removed, while event to be identified 

as a signal must satisfy the dE/dx requirement defined below: 
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• dE/dx requirement: 

The proton must have a measured dE/dx < 9 keV/ cm (see Figure 2). For the accompanied 
e or /i, the measured dE/dx must agree with its expected value for the particle assignment 
and the measured momentum. 

No events in the data sample survived the selection criteria, while 1.2 background events 
were expected from the Standard Model background processes. 



6 Efficiency and systematic errors 

To estimate the detection efficiency of the signal and to investigate possible bias in the selection 
of events, preselection and particle identification criteria were applied to reconstructed Monte 
Carlo signal events. The shaded histogram in Figure [I], with arbitrary normalization, shows the 
distribution of the expected total calorimeter energy, E tot , associated with each of the measured 
tracks of pre-selected Monte Carlo simulated signal events. The low energy peak centered at 3 
GeV comes from the muons, while the peak centered at the beam energy (45 GeV) comes from 



the electrons. The broad peak in the middle of the spectrum comes from the protons [T^, 20]. 
No such central region broad peak is seen either in the data (points) or in the Monte Carlo 
simulated background (open histogram) events as shown in Figure [I]. 

The preselection criteria limited the acceptance of the signal to the barrel region of the 
OPAL detector where the jet chamber provided a large number of hits for each of the pre- 
selected charged tracks for the dE/dx and momentum measurements. These criteria ensured 
a better than 5.3 % resolution in momentum and a 3.5 % resolution on dE/dxfor a 45 GeV 
charged particle. As a result, 3a separation between protons and muons, and 4<r separation 



between protons and electrons at 45 GeV were achieved |10|. The dE/dx distribution of the 
pre-selected data is compared to that of the Monte Carlo simulated background in Figure 0. 
The shaded histogram shows the dE/ dx distribution expected for 45 GeV protons. A dE/dx < 9 
keV/cmcut was used for the identification of protons. 

Table [l] summarizes the event flow history. The equivalent integrated luminosities of the 
Monte Carlo generated background sample sizes are shown in row 1. Row 2 shows the numbers 
of these events, normalized to the recorded luminosity (119pb~ ), surviving the requirements 
of low (charge) multiplicity. In row 3, the numbers of Monte Carlo simulated events and data 
events surviving the preselection are compared. The agreement here is excellent. Not shown 
in the Table are the corresponding numbers for the standard model Monte Carlo simulated 
multi-hadronic decays of the Z°, because we found that none of the 1.3 million Monte Carlo 
simulated Z° — > qq events survived the preselection criteria. In addition, six million Standard 
Model Z° decays were generated with JETSET [|TjJ] and subjected to a preselection at the four- 
vector level. We found no event other than the Standard Model e + e~, fi + ^~ or t + t~ pairs that 
survived the preselection criteria. In Table 1, row 4 shows that: 

the pre-selected data were dominated by (97.4 % ) e + e~ , fi + n~ pairs. 

after preselection, the efficiency for the identification of e + e~ was found to be 99.9 %, and 
94.5 % for 
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after the identified e + e~ and pairs were removed, none of the events in the data sam- 

ple survived the dE/dx requirement selection, while 1.2 events were expected from the 
simulated background. 

the detection efficiencies of the signal events for Z° — > p/i and pe, were found to be: 
e m = 0.353 ± 0.006(stat) and e pe = 0.338 ± 0.007(stat) respectively. 

More details on particle identification for the pre-selected events are shown in Table 0. 
The numbers of tracks surviving the preselection criteria are shown in row 1. The numbers of 
electron and muon candidates identified in the data sample are compared to those of the Monte 
Carlo simulated backgrounds in rows 2 and 3. Proton candidates shown in row 4 are due to 
the mis-identification of (~ 0.5 % ) the pre-selected charged leptons before imposing the final 
dE/dx requirement. This was anticipated because: 



f a high energy muon or electron can sometimes deposit a great deal of its energy inside 
the calorimeters due to bremsstrahlung, e + e~ pair production, or deep inelastic nuclear 
interaction (DIS). These processes are characterized by large uncertainties, small cross 
sections, hard spectra, and large energy fluctuations in the generation of electromagnetic 
and hadronic showers, 

f a high energy lepton can undergo DIS in the ECAL, where only the electromagnetic com- 
ponent of the shower is measured by the lead glass, resulting in Eecal < 35 GeV, 

f hadrons are emitted in r decays. 



These effects, in particular the DIS of the e's and /i's, made the proton identification, with 
E tot cuts alone, difficult. However, after removing the identified e + e~ and pairs, the final 

dE/dx requirement efficiently eliminated every one of the 1249 identified proton candidates in 
the data. The 1.2 background events, all identified as Z° — > p+/i, consisted of 0.5 events coming 
from fi + fi~ and 0.7 events from t + t~ . The background events can be understood as DIS 

of one of the high energy muons, while the t + t~ background events are due to mis- identification 
of Z° — > tit 2 , where T\ —* h and r 2 — > fi. 

Allowing the dE/dx and E to t cuts to vary, based on the analysis of the Monte Carlo simulated 
events, we estimated a 0.71 % systematic uncertainty for the detection of the signal. 

Final-state radiation (FSR) is expected to be different for muons and electrons or protons. 
Changing the particle code from /x to e in the final-state of Z° — ► events (see section 3) 

may have caused an underestimation of the FSR correction on the electrons in the Z° — > pe 
events resulting in a larger than expected acceptance of the signal. Based on a study of the 
Z° — > e + e~ events, we found that 0.47 % of the electrons could have a measured E EC al less 
than 35 GeV. This value was quoted as the systematic uncertainty in the estimation of the 
detection efficiency due to FSR. 

Another source of systematic error came from the 0.41 % uncertainty in the collected lumi- 
nosity |]l9fl. In summary, a quadratic sum of all of the above mentioned errors yielded a 0.94 % 



systematic error on the detection efficiency, e pe , and 0.81 % on e P/Lt . 
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7 Results and conclusion. 



Using a data sample of 5.2 million Z° 
the following upper limits: 



nz 



decays and with candidate event found, we obtained 
■ pe) = 4.6 keV 



and 



r(Z° -> P fi) =4.4keV 



at 95 % confidence level [21 



This is the first time that limits have been obtained on the partial widths of Z° decays 
for processes which simultaneously violate baryon-, lepton-, and fermion- number conservation. 
The above limits are comparable to the existing limits on lepton family number violating Z° 
decays M. 
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Data 
Collected in 
1990-1994 


Monte Carlo 


Background 


signal 


(e+e") 




(T+T-) 


total 


(P,e) 




1. 


Sample size 
used in 
analysis 


119 

(pb v ) 


252 
(Pb" 1 ) 


463 
(Pb' 1 ) 


345 
(Pb" 1 ) 




10,000 

(events) 
100 % 


10,000 

(events) 
100 % 


2. 


Multiplicity 

and Status 
requirements 


1,896,853 


440,027 


148,053 


141,117 


729,197 


92 % 


93 % 


3. 


preselection 


139,328 


73,855 


65,436 


149 


139,440 


48 % 


50 % 


4. 


Particle 
identification 


135,679 
( identified 


73,797 

a 


61,807 

S i 


84 

5° -> e+e" 


135,688 

< 


13 % 


14% 
i + /i~ ) 




( 




identifi 


0.5 

ed 


0.7 

as 


1.2 


33.8 % 
gnal 


35.3 % 



Table 1: Comparison of event flows for data and Monte Carlo simulated events. Numbers 
shown in rows 2, 3 and 4 represent events surviving the selection criteria. The Monte 
Carlo background event numbers were normalized to the luminosity of the data. 
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Data 
Collected in 
1990-1994 


Monte Carlo 


Background 


(e+e") 


(^-) 


(T+T-) 


total 


preselection 


278,656 


147,710 


130,872 


298 


278,880 


Electron 
candidates 


147,527 


147,467 


93 


151 


147,711 


Muon 
candidates 


124,996 


3 


125,142 


41 


125,186 


Proton 
candidates 


1 1 A CI 

i,z4y 


53 


914 


78 


1,045 


Proton 
identified 








0.5 

as 


0.7 

Z ^ 


1.2 
/i + p 



Table 2: Comparison of track flows for data and Monte Carlo simulated events during particle 
identification. Numbers shown in the Table give the number of tracks surviving the 
selection criteria. 
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Figure 1: The measured calorimeter energy (data points), E tot , associated with each well mea- 
sured charged track satisfying the event preselection criteria for the data. The 
Monte Carlo simulated fi + fi~ ,e + e~ , and t + t~ events are shown by the open his- 
togram for comparison. The shaded histogram represents the Monte Carlo simulated 
Z° — > pe/p/i decays with arbitrary normalization. 
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Figure 2: The distribution of specific ionization energy loss, dE/dx, of the well measured 
charged tracks (data points) that survived the event preselection. The open his- 
togram shows the corresponding distribution for the Standard Model Monte Carlo 
simulated background charged lepton pairs. The dE/dx distribution of a 45 GeV 
proton which survived the same selection criteria is shown as the shaded histogram 
with arbitrary normalization. 
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